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Introduction 


1  . 


This  report  describes  work  performed  by  Northwest  Research 
Associates  on  contract  with  the  Air  Force  Geophysics  Laboratory 
( AFGL ) ,  Ionospheric  Effects  Branch  (LIS),  to  provide  hardware  and 
software  support  enabling  the  presentation  of  a  special  scaled 
data  set  characterizing  ionospheric  effects  at  Thule,  Greenland. 

Hardware  support  consisted  of  the  responsibility  of  insuring 
operation  and  maintenance  of  the  monitoring  equipment;  repair  of 
associated  equipment  including  computers,  tape  drives,  receivers 
and  other  components;  implementation  of  system  upgrades;  and 
installation  of  additional  systems. 

Software  support  of  the  TRN  objective  included  installing 
and  customizing  operating  system  upgrades,  installing  and  learn¬ 
ing  new  vendor  software  to  assist  others  on  the  team,  writing 
utilities  and  software  tools,  and  upgrading  existing  routines. 

Additional  activities  included  the  organization  and  cata¬ 
loging  of  tapes  and  charts  coming  from  the  study  site  and  acqui¬ 
sition  of  necessary  materials  required  for  the  study. 

The  study  consisted  of  a  five  step  process: 

(  1  )  Collecting  five  channels  of  raw  data  determined  from 
signals  received  from  GPS  system  satellites  onto  tape  and  chart 
records . 

(2)  Cataloging  these  items  as  they  were  received  at  the 

AFGL. 

(3)  Identifying  elements  of  the  data  for  potential  analyti¬ 
cal  significance  to  the  study. 

(4)  Converting  the  raw  data  and  performing  first-phase 
analysis . 

(5)  Performing  final  analysis  on  data  passing  the  first 
phase.  The  results  of  this  process  were  plotted  and  provided  to 
the  sponsor,  as  required  by  the  TRN. 

2.  Objectives 

The  objectives  of  this  TRN  were  to  provide  engineering  and 
technical  support  for  the  enhancement  and  operation  of  ground- 
based  satellite  monitoring  equipment  that  provide  data  important 
in  the  evaluation  of  ionospheric  effects  on  transionospher ic  RF 
systems . 
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These  entailed: 


a.  Development  and  enhancement  of  the  GPS  data  acquisition 
system  at  Thule  AB,  Greenland,  and  AFGL  at  Hanscom  AFB, 
Massachusetts . 

b.  Installation,  maintenance,  and  operation  of  the  afore¬ 
mentioned  systems. 

c.  Organization,  categorization,  and  analysis  of  the  vari¬ 
ous  information  measures  obtained  from  these  systems. 

d.  Improvement  and  development  of  new  systems  which  facil¬ 
itate  analysis  and  which  will  accommodate  various 
resources  available  at  AFGL  (such  as:  DEC  pdp  11/03's, 
DEC  Micro  pdp  11,  DEC  VAX,  CDC  Cyber,  and  various  PC 
systems ) . 

e.  Utilization  of  existing  systems  and  enhancment  and  cre¬ 
ation  of  new  analysis  systems  to  develop  parameter 
statistics  and  plots. 

3.  Operation  and  Maintenance  of  GPS  Monitoring  Equipment 

One  of  the  key  elements  in  the  GPS  monitor  system,  and  a 
primary  focus  of  attention,  is  the  GPS  satellite  receiver. 
Installation,  calibration,  repair,  and  maintenance  of  the 
receiver  systems  at  Thule  require  one-  to  two-week  trips  every 
two  to  three  months. 

Often  it  is  necessary  to  ship  equipment  back  to  the  labora¬ 
tory  at  AFGL  for  repair.  Other  situations  are  dealt  with  by 
phone,  analyzing  the  problems  and  instructing  the  Danish  techni¬ 
cians  how  to  proceed. 

In  November,  1987,  an  on-site  receiver  evaluation  at  Thule 
revealed  the  necessity  of  replacing  a  receiver  component  for 
which  there  was  no  readily  available  spare.  While  there,  a  new 
receiver  arrived  at  AFGL,  prompting  the  decision  to  curtail  the 
trip,  bench-check  the  new  receiver  at  AFGL,  and  obtain  the  neces¬ 
sary  replacement  component.  Another  trip  to  Thule  was  made  to 
install  the  new  receiver  and  repair  and  retain  the  old  receiver 
on  site  as  a  backup. 

In  early  1  988  a  decision  was  made  to  expand  observation  at 
Thule.  A  second  system  was  assembled  from  spare  parts  which 
would  allow  observation  of  high-elevation  satellites  on  the  cur¬ 
rent  system  and  concurrent  observation  of  low-elevation  satel¬ 
lites  on  the  newly  assembled  system.  In  early  June  the  new  sys¬ 
tem  was  installed,  and  upgraded  acquisition  software  was 
installed  on  both  systems  at  that  time. 

The  GPS  development  system  at  AFGL  was  deployed  to  Kwajalein 
Island  in  the  South  Pacific  in  late  July,  1988,  to  participate  in 
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the  DNA-sponsored  Propagation  Effects  Assessment-Kwa jalein  (PEAK) 
campaign.  NWRA  personnel  installed  and  operated  the  system  for 
the  six-week  duration  under  TRN  #13.  The  system  was  subsequently 
shipped  to  Thule  to  participate  in  another  campaign  to  be  con¬ 
ducted  in  December,  1988,  and  to  serve  as  a  complete  backup  for 
the  two  systems  on-site. 

4.  Software  Upgrades  and  Enhancements 

A  preparatory  effort  was  undertaken  at  the  outset  of  the  TRN 
to  implement  refinements  and  anticipate  potential  problems.  This 
entailed  upgrading  and  enhancing  the  software  portions  of  the  GPS 
collection  and  analysis  systems,  and  correcting  and  improving  the 
hardware  portions  of  the  GPS  acquisition  system.  Additionally, 
software  tools  were  created  to  extend  ability  to  deal  with  known 
data  requirements. 

4.1  Data  Acquisition  Software 

The  data  acquisition  system  software,  which  operates  on  the 
DEC  pdp  11/03's,  was  analyzed  and  upgraded.  The  satellite 
'window'  portions  were  streamlined,  and  enhancements  were  devised 
to : 

a.  Improve  satellite  'window'  information  entry,  deletion, 
and  revision. 

b.  Provide  for  automatic  update  of  satellite  'window'  ele¬ 
ments  that  previously  had  become  gradually  obsolete 
over  a  few  weeks.  This  eliminated  a  cause  for  frequent 
attention  to  satellite  'window'  elements. 

c.  Offer  a  choice  of  Danish  instructions  as  an  alternative 
to  English  to  reduce  the  possibility  of  misunderstand¬ 
ing  by  the  technicians  who  perform  day-to-day  opera¬ 
tions.  The  translations  were  performed  by  NWRA's  Mr. 
C.  Charley  Andreasen. 

4.2  Improvements  to  the  Real-Time  Routine 

a.  The  real-time  routine  that  acquires  data  was  stream¬ 
lined  as  well.  The  following  corrections  were  made: 

( 1  )  Solving  an  intermittent  failure  to  detect  the  end 
of  a  satellite  pass,  which  caused  the  system  to 
lose  useful  data  over  the  succeeding  24  hours. 

(2)  Performing  a  check  to  prevent  an  operator  from 
erroneously  instructing  the  system  immediately  to 
start  taking  data  from  a  satellite  that  is  not 
available  at  that  time. 
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(3)  Allowing  the  system  to  detect  the  loss  of  'lock* 
by  the  receiver  on  the  current  satellite,  thereby 
halting  acquisition  of  invalid  data. 

b.  Enhancements  to  the  real-time  routine  included: 

( 1 )  Creation  of  a  Danish  language  version  of  the 
acquisition  routine,  automatically  invoked  accord¬ 
ing  to  the  choice  made  in  the  'window'  routines 
(Separate  English  and  Danish  version  of  this  rou¬ 
tine  were  necessary  due  to  its  size  and  hardware 
restrictions ) . 

(2)  Implementation  of  a  header  that  records  the 
satellite  (space  vehicle)  number,  Julian  date  and 
year  of  the  pass,  and  the  session  operator. 
(Previously,  the  satellite  number  and  Julian 
date/year  were  determined  and  added  manually  dur¬ 
ing  the  analysis  phase.)  Automation  eliminated 
the  potential  for  error  and  reduced  manpower 
requirements.  Recording  the  session  operator  pro¬ 
vided  a  means  for  tracking  operator  errors. 

4.3  Problem  Detection  and  Recovery 

Software  tools  were  created  to  facilitate  problem  detection 
and  recovery,  which  included: 

a.  A  routine  to  allow  a  quick  look  at  raw  data,  averaged 
to  six-second  quantities  (data  are  collected  at  20  val¬ 
ues  per  second  in  binary  format),  enabling  trace-back 
of  faults  noted. 

b.  A  routine  to  select  valid  data  from  a  satellite  pass 
file  that  failed  to  end  appropriately  and  save  it  in  a 
separate  file. 

c.  A  program  to  copy  selected  files  to  a  separate  tape  to 
enable  analysis  of  only  those  files. 

4.4  Primary  Analysis  Routine 

Refinements  were  made  to  the  primary  analysis  routine,  which 
converts  the  raw  GPS  measurements  into  slant  Total  Electron  Con¬ 
tent  (TEC)  values,  enhancing  functionality.  Access  to  the  header 
information,  now  recorded  in  the  acquisition  system,  was  added. 
Multipath  correction  of  Differential  Group  Delay  procedures, 
developed  under  TRN  #13,  was  added.  Other  enhancements  were 
attempted,  but  architectural/system  constraints  of  the  CDC  Cyber, 
the  primary  analysis  system,  prevented  their  implementation. 
Alternative  solutions  are  under  consideration. 
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4.5  Quick-Look  Summaries 


Procedures  were  established  for  the  generation  of  'quick- 
look'  data  summaries.  An  activity  quantifier  measure  (see  the 
Table  on  page  6)  was  established  whereby  hourly  activity  levels 
in  GPS  TEC  measurements  and  polar  beacon  satellite  250-MHz 
amplitude  scintillation  measurements  were  tabulated.  A  visual 
examination  of  chart  records  from  Thule  would  be  performed  as 
they  arrived  and  a  tabulation  made  accoraing  to  the  standard 
developed.  This  enabled  analysis  to  be  focused  on  pass  files 
with  significant  ionospheric  activity. 

5.  Production 

Production  work  and  preparatory  work  were  performed  concur¬ 
rently  in  the  early  stages  of  the  TRN.  Production  work  consisted 
of  reviewing  the  tapes  and  charts  arriving  from  Thule,  cataloging 
them,  visually  examining  the  chart  records,  and  tabulating  the 
activity  levels  according  to  the  criteria  described  above.  Using 
criteria  outlined  in  the  Table,  active  periods  of  significant  TEC 
gradients  or  significant  fades  from  250-MHz  beacon  satellite  mea¬ 
surements  were  identified.  Raw  data  occurring  during  these 
active  periods  were  then  processed  by  the  first  level  of  GPS 
analysis,  which  converted  GPS  differential  carrier  phase  measure¬ 
ments  to  slant  TEC  variation. 

Differential  carrier  phase  measurements  were  considered  the 
best  available  measure  to  indicate  the  rapidity  of  change  in  the 
ionosphere.  Absolute  TEC  measurements,  available  from  a  measure 
of  differential  group  delay  output  from  the  GPS  receiver,  are 
less  suitable  for  this  study  because  of  susceptibility  to  system 
noise  and  multi-path  effects. 

To  perform  the  final  analysis,  an  'event'  analyzer  was 
developed.  This  enabled  large  quantities  of  'significant'  activ¬ 
ity  to  be  efficiently  analyzed  and  quantified  according  to 
activity  level. 

Several  approaches  were  examined  and  tested  for  characteriz¬ 
ing  TEC  measurements  specific  to  a  particular  satellite  pass.  It 
was  decided  to  assess  the  TEC  data  by  determining  a  measure  of 
the  TEC  gradient  over  a  sliding  360-second  sampling  period.  The 
measured  'event'  gradient  was  determined  by  the  difference  of  the 
means  of  five  samples  af  40-second  intervals  in  the  first  half 
and  the  second  half  of  the  sampling  period.  An  'event'  level  of 
'1'  was  defined  as  a  measured  difference  of  3.3x1 0^  electrons/m2 
occurring  in  the  40-second  time  difference  interval. 

A  routine  was  developed  to  tabulate  slant  TEC  measurements, 
'event'  magnitude,  and  satellite  azimuth  and  elevation  over  the 
duration  of  the  pass.  These  routines  were  further  upgraded  to 
determine  the  probability  of  an  event  occurrence  and  the  proba¬ 
bility  of  no  event  during  a  period  twice  the  defined  'event' 
length  and  display  the  results.  Twelve  representative  samples  of 
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Table  1.  Activity  quantifier  for  GPS  and  250  MHz  beacon  data. 


GPS 

250  MHz 
Beacon 

Comment 

Quiet 

1 

A>0.5  TEC  u/ 1 0  min 

SI > 2  dB 

No  detectable  activity. 

Small 

2 

A<  0 . 5  TEC  u/5  min 

2<S I < £  dB 

Detectable  TEC  variation, 
detectable  scintillation. 

Moderate 

3 

A>1  TEC  u/5  min 

5<_SI<9  dB 

Smooth  rolling  structures. 

Disturbed 

4 

A< 1  TEC  u/30  sec 

1  0_^S  I  <  1  5  dB 

Sharp,  short  gradients. 

Very 

Disturbed 

5 

A< 1  TEC  u/30  sec 

SIM  5  dB 

Very  sharp,  sustained 
gradients . 

the  data,  averaging  about  2  samples  per  month,  are  presented  in 
Appendix  A,  covering  the  analysis  period. 

Amplitude  scintillation  data  at  250  MHz  was  categorized  and 
tabulated  and  compared  wxth  simultaneous  TEC  data  provided  by  the 
'event'  analyzer.  A  standard  procedure  was  applied  to  derive 
fade  depth  measures  from  the  charts  at  fifteen-minute  intervals. 
The  procedure,  presented  in  Appendix  B,  was  identical  to  that 
used  by  AFGL  to  qenerate  fade  depth  occurrence  statistics  since 
1979. 


A  database  structure  was  constructed  to  allow  accumulation 
of  event-occurrence  data  according  to  various  parameters.  The 
data  selected  from  the  'event*  analyzer  were  incorporated  into 
it,  as  was  the  probability  of  non-occurrence  of  events  as 
described  above.  Additional  data  include  tabulated  250  MHz  fade 
depth  values,  solar  and  magnetic  indices,  and  several  other 
parameters . 

6.  Final  Data  Analysis 

The  database  utilized  for  final  analysis  included  approx¬ 
imately  130  records  (pass  files)  from  October  1987  through  March 
1988.  After  culling  of  incomplete  records,  about  100  records 
were  available  for  analysis. 

The  AFGL  measurement  data  base  obtained  during  the  reporting 
period  is  summarized  in  Appendix  C  and  is  entitled  "GPS  Sunspot 
Activity  Special  Report."  Tabulated  data  are  presented  that 
show,  for  corresponding  record  ID's,  'event'  levels  as  a  function 
of  observed  fade  depth.  As  shown,  an  additional  field  was  added 
to  the  database  called  'QA'  (quality  assurance),  which  denoted 
whether  the  record  was  valid  (V)  for  analysis  or  invalid  (I). 

The  records  were  retained,  nonetheless,  to  allow  correction  of 
their  deficiencies  for  inclusion  at  a  later  time. 

Final  analysis  involved  the  statistical  evaluation  of  valid 
database  entries.  This  required  the  development  of  a  set  of  sta¬ 
tistical  analysis  programs  written  in  DBASE  III  programming  lan¬ 
guage  for  use  on  the  Zenith  PC. 

The  following  data  presentations  were  produced: 

a.  Curves  of  the  probability  of  a  period  without  an  event 

for  all  active  data  in  the  dataset  are  presented  in 

Appendix  D.  Curves  of  the  probability  of  a  period 

without  an  event  are  shown  for  each  discrete  fade  depth 
category  observed.  The  first  part  of  the  appendix  pre¬ 
sents  the  fade  depth  algorithm  and  indicates  several 
methods  of  categorizing  the  data. 

b.  Histograms  of  the  number  of  hours  and  percentage  of 
occurrence  of  fade  depth  levels  within  active  TEC  and 
all  TEC  5l  ervations  are  presented  in  Figure  1  and 
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through  March  1988. 
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Figure  2  respectively.  Active  TEC  corresponded  to 
those  times  when  quick-look  analysis  indicated  that  the 
'event'  level  was  3  or  greater. 

c.  Histograms  combining  monthly  summaries  of  diurnal 
occurrence  of  TEC  activity  and  250-MHz  amplitude  scin¬ 
tillation  over  the  study  period  are  presented  in  Figure 
3.  Review  of  these  data  showed  that  times  of  TEC 
activity  corresponded  to  times  of  active  250-MHz  ampli¬ 
tude  scintillation.  A  specific  case  study,  summarized 
in  Figure  4,  highlighted  a  fortuitous  intersection  of 
the  GPS  and  polar  beacon  look  directions  during  high 
TEC  activity.  This  study  showed  that  there  was  a 
direct  relationship  between  the  observed  TEC  enhance¬ 
ment  and  250-MHz  scintillation. 

d.  Histograms  are  presented  in  Figure  5  that  show  the 
occurrence  of  15-minute  fade  depth  indices  from  1979 
through  1984.  Data  are  broken  down  into  categories 
comparable  to  Total ,  Active ,  and  Exceeding  Active . 

Results  from  the  histograms  of  fade  depth  occurrence  during 
the  study  periods  showed  that: 

a.  The  fade  depth  activity  levels  occurring  throughout  the 
entire  reporting  period  corresponded  to  the  very  lowest 
monthly  activity  levels  that  occurred  during  the  previ¬ 
ous  (1979-1980)  solar  maximum,  and  are  tabulated  in 
Figure  5. 

b.  These  data  indicate  that  the  fade  depth  activity  levels 
occurring  during  periods  of  active  TEC  of  the  study 
corresponded  to  less  than  those  during  the  lowest 
activity  month  observed  during  the  previous  solar  maxi¬ 
mum. 

c.  A  study  was  made  comparing  the  occurrence  of  15-minute, 

250-MHz  fade  depth  in  the  data  set  obtained  during  the 
current  study  period  with  published  statistics  of  fade 
occurrence  over  the  period  1979  through  1984.  Fade 
depths  were  manually  determined  from  chart  records  of 
250-MHz  intensity  scintillation  using  the  identical 
process  applied  to  obtain  published  fade  depth  statis¬ 
tics  (see  Appendix  B).  As  shown  in  Figure  5,  the 
greatest  intensity  fade  levels  occurred  during  the  pre¬ 
vious  sunspot  maximum.  From  a  review  of  the  data  it 
can  be  seen  that  in  any  given  month  the  second  highest 
fade  and  highest  fade  category  can  occur  up  to  60  per¬ 
cent  and  45  percent  of  the  time  respectively.  The 

occurrence  of  fade  depth  during  the  current  study 
period  summarized  in  Figure  1  reveals  that  the  highest 
fade  depth  was  not  observed  to  occur. 
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Fade  DATA  QUANTITY 

Percent 


Figure  2.  Percent  occurrence  of  various  intensity  fade  levels 
within  all  TEC  observations,  from  Thule,  Greenland, 
October  1987  through  March  1988. 
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Figure  3.  Monthly  summaries  over  the  study  period  showing 
occurrence  of  TEC,  Scintillation  Activity,  and  foF2  as 
a  function  of  universal  time. 


foF2  foF2  foF2  foF2 


Figure  4.  Measurements  from  GPS  and  polar  beacon  satellite 
passing  within  1  degree  show  TEC  variations  linked  to 
strong  amplitude  scintillations  at  250  MHz. 


12 


3 


o 


o 


o 


3  3  3  3 


siaaai 


*a  rt  o  2 


Ik  At  M 


03 

tA 

c0 

p» 

P> 

03 

00 

iO 

t/J 

CO 

X> 

— 

co 

00 

P> 

J- 

00 

CO 

CO 

■o 

s: 

u 

00 

i 

CO 

CO 

as: 

— 

fNk 

00 

i 

J- 

• 

cO 

r^ 

2e: 

«p> 

00 

o 

— 

i 

o 

4-> 

E 

U 

co 

2 

o 

O- 

E 

4-> 

o 

U 

■P 

t- 

O 

0) 

c+- 

<A 

E 

cO 

o 

CO 

4-> 

v 

P 

cO 

«0 

-o 

"O 

C0 

tu 

P> 

-J 

:> 

cO 

i— » 

T3 

h- 

H 

O 

O 

Ui 

H 

c 

> 

HM 

o 

o 

H 

4J 

p 

O 

<C 

03 

03 

r* 

cr 

3 

X> 

c 

CO 

CO 

t» 

U 

XI 

CO 

cO 

0) 

a. 

§ 

1 

03 

U 

X 

o 

u 

0) 

cO 

CO 

c0 

P 

P 

P 

CO 

CO 

CO 

O 

o 

o 

X 

<0 

e 

L. 

•O 

O 

tA 

4-> 

<o 

s. 

<TJ 

<U 

>» 


X  <A 
fO  -C 

E  -M 

C 

u  o 
<0  E 

o  vo 

to 

e» 

4->  U 
CO  C 
0) 
i-  s- 

03  t- 
<U  =3 

>»  o 
u 

<4-  O 

o 

4J 
(O  C 

jc  <u 

4->  U 
C  V- 
Q  QJ 


00  "O 

c 
i-  co 
o 

V-  >3 

p> 
to  T* 

xj  > 

<U  *r» 

a>  40 
u  u 
x  co 

c 

U  *r- 
O 

to 
to  T3 
r-  01 
10  03 
3  U 

cr  x 


>◄ 


2  2 


33M3V«n330  39V1N33M34 


Appendix  A: 


Twelve  representative  examples  of  equivalent  slant  TEC 
covering  period  day  290,  1987,  through  day  88,  1988 
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Appendix  B: 


Description  of  the  standard  AFGL  procedure  to  measure  fade 
depth. 
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Instructions  for  Measuring  Beacon  Signal  Strength  Variation 

Beacon  signal  strength  can  measure  ionospheric  fading  of  the  satellite  signal  during 
ionospheric  disturbances.  Beacon  signal  strength  is  recorded  concomitantly  with  the  several 
measurements  used  by  GPS  to  determine  Total  Electron  Content  (TEC)  and  written  to  a 
chart  record. 

To  determine  a  measure  of  beacon  signal  strength  variation,  a  measuring  device 
must  be  made.  A  step  calibration  of  the  signal  is  normally  performed  daily  by  the  techni¬ 
cians  at  the  receiver  site.  The  measuring  device  can  be  made  from  the  chart  record  of  this 
calibration.  Attached  is  a  sample  of  such  a  calibration  (attachment  1).  It  will  be  referred  to 
in  the  following  explanation  of  how  to  create  the  measuring  device. 

Although,  the  step  calibration  of  the  signal  strength  is  usually  done  once  per  day,  this 
is  not  always  possible.  Therefore  a  step  calibration  which  preceeds  the  data  to  be  analyzed 
should  be  used  to  make  the  measuring  device  rather  than  one  which  succeeds  the  data  even 
though  the  succeeding  step  calibration  may  be  closer  in  time  to  the  data  under  considera¬ 
tion.  (This  is  because  the  subsequent  calibration  may  effect  altered  receiver  tuning.) 

A  3x5  card  is  the  easiest  item  to  use  as  a  measuring  tool.  Place  the  top  edge  of  the 
card  along  the  top  graph  line,  of  the  signal  strength  channel,  which  runs  above  the  calibra¬ 
tion  marks  (see  A  in  attachment  1).  Place  it  so  the  right  edge  abuts  the  first/highest  step  of 
the  stair-step  calibration  marks  (see  B).  Make  a  line  on  the  right  edge  of  the  card  which 
corresponds  to  that  highest  mark.  Identify  that  mark  as  'O'.  Now,  advance  the  card  along 
the  top  line  (A)  to  the  second  step  from  the  first  mark  (C).  Again,  make  a  mark  on  the  card 
and  label  this  '2'.  Repeat  this  for  '4'  and  '6'  (see  D  &  E).  After  these,  slide  the  card  until  its 
edge  covers  the  long  pulse  step  (see  F).  Make  a  mark  on  the  card  which  corresponds  to  a 
point  halfway  down  from  that  step  and  the  next  and  label  that  '8'. 

The  calibration  changes  after  this.  Marks  are  made  at  every  step  rather  than  at 
alternate  steps.  Each  of  these  steps  represents  an  increment  of  two.  For  example,  the  card 
is  now  advanced  to  the  next  step  (see  G),  marked  and  labeled  10'.  The  position  after  that 
(see  H)  is  labeled  12'.  Repeat  this  until  '24'  (see  I). 

The  creation  of  the  measuring  device  is  now  complete.  The  zero  mark  is  equivalent 
to  66dB;  the  twenty-four  mark  is  equivalent  to  a  range  of  about  91dB. 

Locating  a  desired  pass  file  for  analysis  can  vary  in  degree  of  effort  depending  on  the 
annotation  of  the  corresponding  chart  record.  Minimal  annotation  usually  consists  of  a  time 
and  date  written  at  the  24  hour  delineation  and/or  a  time  written  at  the  step  calibration. 

Time  is  denoted  on  a  chart  record  by  a  continuous  straight  ink  line  with  tick  marks  at 
regularly  spaced  intervals.  The  tick  marks  vary  in  type;  there  are  three  (see  attachment  2). 
One  type  is  the  minute  mark  (see  A);  another  type  is  the  ten-minute  mark  (see  B);  the  last 
mark  is  the  hour  mark  (see  C). 
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Assuming  the  chart  record  is  minimally  annotated,  to  locate  the  desired  GPS  pass  file 
flip  through  the  chart  record  page  by  page  until  a  24  hour  annotation  or  step  calibration  that 
immediately  preceeds  the  desired  period  is  found.  Now  locate  the  hour  mark  after  that. 
Count  the  hours  until  the  one  preceeding  the  pass  file  start  is  reached.  (Hour  marks  are 
approximately  in  the  same  place  on  alternate  pages.)  Between  this  hour  mark  and  the  next 
a  GPS  calibration  mark  in  the  chart  record  will  be  apparent  on  one  or  several  channels  of 
the  chart  record  (see  D).  This  cal  mark  denotes  the  beginning  or  end  of  a  GPS  pass  file. 
From  this  mark  until  the  succeeding  one  is  the  extent  of  a  particular  GPS  pass  file. 

Beacon  signal  strength  is  determined  by  locating  the  largest  events  on  the  Beacon 
channel  in  fifteen  minute  intervals  and  measuring  their  widths.  Define  fifteen  minute  inter¬ 
vals  from  the  first  ten  minute  tick.  An  interval  should  begin  at  every  alternate  ten  minute 
tick.  Observe  the  number  of  ticks  from  the  first  GPS  cal  mark  to  the  first  ten  minute  tick.  If 
the  number  of  minutes  is  seven  or  greater,  determine  a  beacon  signal  strength  for  this  inter¬ 
val  (include  the  area  prior  to  the  cal  mark  so  fifteen  minutes  are  used).  Determine  a  signal 
strength  value  for  each  interval  until  the  ending  cal  mark  is  reached  or  passed. 

Measuring  the  signal  strength  is  accomplished  by  the  following  method.  Observe  a 
fifteen  minute  interval.  Find  the  third  highest  signal  level  in  the  interval  and  draw  a  line  at 
its  peak.  Next,  locate  the  third  lowest  signal  level  in  the  interval  and  draw  a  line  at  its  nadir. 
Take  the  3x5  card/measuring  device  created  earlier,  align  its  edge  with  the  top  graph  line  of 
the  channel,  near  the  mark  made  for  highest  signal  level.  Note  the  corresponding  value  on 
the  card.  Repeat  this  for  the  nadir  mark.  Again,  note  the  corresponding  value.  Calculate 
the  difference  between  these  two  values  and  note  it  on  the  chart  record.  This  is  the  value 
fo.  Beacon  signal  strength. 

When  the  positive  excursions  of  the  Beacon  signal  strength  exceed  both  the  'zero' 
mark  on  the  measuring  card  and  the  top  graph  line  (A),  one  additional  dB  unit  should  be 
added  to  the  Beacon  signal  strength  value  determined.  When  several  negative  excursions 
reach  the  'twenty-four'  mark  on  the  measuring  card  (minimum  that  the  calibration  can 
qualify)  then  2dB  should  be  added  to  the  Beacon  signal  strength  value.  When  'many'  nega¬ 
tive  excursions  reach  the  'twenty-four'  mark,  then  four  or  more  dB  should  be  added. 
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Beacon  Scintillation  Data  Analysis 
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Appendix  C: 

Database  utilized  in  final  data  analysis. 
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Page  No. 
09/14/88 
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GPS  SUNSPOT  ACTIVITY 
Special  Report 


REC 

JULIAN 

UT 

UT 

10 

OAY  DATE 

START 

STOP 

1 

281  10/08/87 

171754 

203434 

2 

282  10/09/87 

171309 

203036 

3 

282  10/09/87 

203343 

220050 

4 

282  10/09/87 

221134 

233135 

S 

286  10/13/87 

1  70S  l  7 

225921 

S 

290  10/17/87 

124446 

151048 

7 

290  10/17/87 

151500 

165033 

8 

292  10/19/87 

165002 

200643 

9 

293  10/20/87 

162110 

1  75841 

10 

293  10/20/87 

1 80205 

19383S 

1 1 

293  10/20/87 
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223937 

12 

294  10/21/87 

14442 

2S857 

14 
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62704 

903SI 

IS 
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134847 

16 
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1721 17 

1 85848 

17 
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18 

307  11/03/87 
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214933 

19 

307  11/03/87 

225243 

3842 

20 

308  11/04/87 

4212 

15842 

21 

309  11/05/87 

195849 

213750 

22 

309  l l /OS/87 

224640 

2641 

23 
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Appendix  0: 


Curves  of  the  probability  of  a  period  without  an  event  for  ail 
active  data  given  several  methods  of  categorizing  the  data. 


Peak-to-Peak  Values  Measured 


0  _> 
8.6 

15.7 

22.8 
28.0 


Method  1  : 


Method  2: 


Method  3: 


Method  4 : 


Method  5 : 


Corresponding 
Fade  Depth  Value 


’  to  P  <  8 . 6 

0 

dB 

P  to  P  15.7 

5 

dB 

’>  P  to  P  <  22.8 

1  0 

dB 

>  P  to  P  <  28.0 

1  5 

dB 

<  P  to  P 

29 

dB 

(Usually  5-12  P  to  P  numbers  per  pass.) 


Fade  Depth  Analysis  Algorithms 


Take  mean  of  all  peak-to-peak  numbers  in  a  pass  and 
determine  corresponding  fade  depth  value  to 
characterize  the  period; 

Determine  corresponding  fade  depth  value  for  each 
peak-to-peak  value  in  a  pass,  sort  fade  depth  values 
in  ascending  order  and  take  value  just  greater  than 
half  as  characterizing  the  period; 

Same  as  Method  2  except  take  value  just  greater  than 
upper  1/3  as  characterizing  the  period: 

Determine  the  highest  peak-to-peak  value  for  the  pass 
period  and  determine  the  corresponding  fade  depth 
value  to  characterize  the  period; 

Determine  corresponding  fade  depth  for  all  peak-to- 
peak  values  and  use  all  of  them  to  characterize  the 
period . 


EVENT  MAGNITUDE  EVENT  MAGNITUDE 
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FADE  DEPTH  LEVEL  ALL 


METHOD  1 
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METHOD  1  METHOD  1 

FADE  DEPTH  LEVEL  5  FADE  DEPTH  LEVEL  10 


NT  MAGNITUDE 


0-0  0.5  1.0  0.0  0.5  1.0 

PROBABILITY  OF  NO  EVENT  PROBABILITY  OF  NO  EVENT 


METHOD  2 

FADE  DEPTH  LEVEL  10 


METHOD  2 

FADE  DEPTH  LEVEL  15 
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EVENT  MAGNITUDE 


METHOD  3 

FADE  DEPTH  LEVEL  0 


METHOD  3 

FADE  DEPTH  LEVEL  5 


METHOD  3 

FADE  DEPTH  LEVEL  10 


METHOD  3 

FADE  DEPTH  LEVEL  15 
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